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Abstract 



Starting from the general cosine form for the axion effective potential, we quantize the axion 
and show that the result is described by a naturally flat potential, if interactions with other 
particles are not considered. This feature therefore restores the would-be Goldstone-boson 
nature of the axion, and we calculate the corresponding vacuum energy density, which does not 
need to be too large by orders of magnitude compared to Dark Energy. 



Introduction. It is known that a quantized scalar theory must be described by a convex one- 
particle-irreducible (IPI) potential [T], as a consequence of the so-called spinodal instability, 
where the restoration force on field fluctuations vanishes. In the situation of a bare potential 
presenting a concave region, as a Mexican-hat potential, the corresponding IPI potential can 
be obtained by a procedure similar to the Maxwell construction in Thermodynamics [2j. This 
property is valid for a self-interacting scalar field only, and therefore does not apply to a situation 
where the scalar also interacts with other degrees of freedom, such as the Higgs boson in the 
Standard Model. 

An example of the self-interacting case is provided by the Sine-Gordon model, in 1+1 dimen- 
sions, where Wilsonian renormalization studies have been done, dealing with the corresponding 
spinodal instability [3]. In these studies, the authors consider the coarse grained potential, 
which interpolates the bare potential and the IPI potential, when the cut off of the theory 
decreases. They find that, between two minima of the cosine bare potential, the coarse grained 
potential flattens when the cut off decreases, and becomes fiat in the IR limit. We note here 
that this flatness holds in a wider area than the one where the bare potential is concave, and 
the flat IPI potential joins the minima of the bare potential. As a consequence, starting from a 
cosine bare potential, the only way for the IPI potential to be convex is to be flat. Other studies 
involving a bare cosine potential, and the corresponding flattening due to spinodal instability, 
have been made in the context of quitessence Pj. The authors consider a description based on 
spinodal decomposition, using a potential function of two variables: the mean scalar field and 
its fluctuation condensate, in the framework of a Hartree approximation. 

We consider here a cosine potential, in 3+1 dimensions, which is characteristic of axions, 
and we show how a flat potential indeed arises. As a consequence, starting with a cosine 
bare potential, the axion self-interactions gradually flatten the potential, and we are left with 
a vacuum energy only, which we calculate. A novel feature of the present study is that we 
quantize the axion degrees of freedom as compared to regarding the axion as a classical field in 
conventional cosine-based expressions for the potential. Using this fully quantized treatment 
we derive an expression for the self-interacting axion vacuum energy and we see that in some 
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heavy axion models, this vacuum energy is not necessarily orders of magnitude larger than 
observed and estimated values for Dark Energy. Coupling the axion to other particles of the 
Standard Model, though, can prevent the flattening of the potential, and we discuss a specific 
example, where the axion is coupled to fermions. We derive conditions to avoid the spinodal 
instability, such that no flattening of the axion potential occurs. 

This article is structured as follows. We first review the quantization of the axion model, 
in order to define our notations, and we show that the corresponding IPI potential must be 
convex. We then derive an exact self consistent equation for this IPI potential, in the form 
of a differential Schwinger-Dyson equation, which is non-perturbative. A fiat IPI potential is 
indeed a solution of this equation, and we calculate the corresponding vacuum energy density. 
Finally, We discuss the coupling of the axion to other particles. 

The axion and its quantization. In the original paper [5], the axion is introduced to build 
a CP conserving theory, from a model with massive fermions coupled to a non-Abelian gauge 
field, where a CP violating term arises at one-loop due to non-perturbative QCD instanton 
effects, [TU]. This original model was refined following experimental and cosmological evidence 
and [7] proposes a more general Lagrangian with three potentially non-zero coupling constants 
from which the original work [5], or invisible axion models such as [11], can be constructed. In 
this paper we use the term axion in this general sense, defined as initially arising as a phase 
- a Goldstone mode - of a complex scalar field, $ which acquires a vacuum expectation value 
($) = /e*^ due to spontaneous breaking of a global U{1) symmetry at a scale /. Explicit 
symmetry breaking occurs at some scale /i, with /i < / at which point the axion can acquire 
mass. Motivated by the need to solve the CP problem, we require that the axion model possess 
a non-zero anomalous gluon coupling term 6GG. References [Zj, [8j, and |9] quote the commonly 
used general expression for the axion potential as K{1 — cos 6), where K has mass dimension 4. 
In arriving at this expression the initial step is integrating over the QCD degrees of freedom. 
The total partition function is 



where represents the QCD degrees of freedom and $ = pe*^. S0 defines then the effective 
action for the axion, which we will quantize. The calculation of Sq is necessarily approximate, 
and an analysis of the form of the axion potential in light of new understandings of the QCD 
vacuum, using supersymmetric gauge and brane theories, is presented in |9j . The authors note 
that the K{1 — cos 6') form for the potential of the QCD axion needs not be limited to one 
cosine dependence, but may have higher harmonics and in general is a smooth periodic function 
of 9 with period 27r. They conclude that with axion models using light quarks mg <^ ^qcd one 
arrives at a cut-off independent potential of the form K{1 — cos 6). 

Motivated by these results and the general form of the Lagrangian presented in [7], we use 
as our general form for the axion action as follows. 



where a„ are coefficients to be determined (units quartic in mass). The action ([2]) can be 
understood as a consequence of the upgrading of the phase of ($) to a dynamical field, that we 
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wish to quantized. An example of quantization of the axion is done in [B], where the authors 
consider quantum fluctuations of the axion, leading to a Bose-Einstein condensate which could 
play the role of cold Dark Matter. 

In terms of the canonically normalized field (j) = f6, this theory is not renormalizable, and 
we will need to keep the cut off A throughout the study, which will be considered a parameter 
of the theory, as well as /. Also, by definition of the cut off, we will always assume that / < A. 
(For completeness, we also assume that Aqcd < A.) 

Quantization of the model (|2]) is based on the partition function. 

Z[j] = e-^[^1 = J V[9] exp (^-Sg - J d^x , (3) 
from which the classical field is defined as 

= ^ = {e) , (4) 

where 



<■■■>= 1 y V[e]{- ■ ■ ) exp (^-Se - j d'x jO^ . 



(5) 



A second functional derivative with respect to the source gives 

6^W 



dclOch (6) 



SjSj 

which is the variance of the set of variables 6{x), therefore positive, such that 

6^W 



< 0, (7) 



which will be used to show the convexity of the IPI potential. The proper graph generator func- 
tional T[6d] is defined as the Legendre transform of the connected graphs generator functional 
W[j]: 

mi] = W[j] - j d^x J 9d, (8) 
where the source j is understood as a functional of 6d, and we have then 



The convexity of the IPI potential U follows directly from the previous properties [Ij. f/ is the 
non derivative part of F: if 6*0 is a constant classical field, F [6^0] = / d'^xU{9o). The corresponding 
source, denoted jo, is also constant, such that the relation between second functional derivative 
in eqs.(l9]) gives then, after integration over space time, 

d^U f . f5^WY^ , , 
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such that, taking the inequahty into account, we have, for any value of 9o, 



and the potential is necessarily convex. As a consequence, if the bare potential is concave, 
at least in a certain range of field values, quantum corrections necessarily contain tree-level 
contributions, in order to compensate the concave features of the bare potential and make the 
IPI potential convex. 

If the axion is coupled to other degrees of freedom, then has to be understood as a matrix, 
with rows and columns corresponding to the different degrees of freedom. The second derivative 
of the IPI potential is then an element of the inverse matrix {6W)~^, which is not necessarily 
positive. In this latter situation, concave features of the bare potential can be stable, and 
quantum corrections are perturbative only. 

Note that, if one considers a periodic axion field, which is restricted to the interval [0; 2it] in 
the path integral ([3]), the latter should be defined with 



but the derivation of the property ([7]) would not change. Indeed, this inequality was not 
obtained by any explicit path integration, but corresponds to a generic property of Z[j]. For the 
same reason, the existence of non-trivial vacua and topological features in the classical theory 
does not alter the property ([7]): all the corresponding information is contained in the path 
integral. Only a saddle point approximation would have to take into account these topological 
features, but no approximation or expansion is made here. 

Finally, by construction, the effective potential must be differentiable, since it is obtained by 
integrating over the original field 6. The only way for the effective potential to be convex, 
differentiable and periodic is to be flat: we will see now that this solution is indeed satisfied. 

Self consistent equation for the IPI potential. Our approach consists in studying the evo- 
lution of the proper graph generator functional F with the amplitude of the decay constant /. 
This functional approach is inspired from the original version given in [12], and will lead us to 
an exact partial differential equation for the IPI potential U. In a first stage, this procedure 
can be seen as a mathematical tool, to obtain a non-perturbative self-consistent equation for 
U. Nevertheless, the physical meaning of the approach is the following. Suppose we start from 
f » m, the interaction is then negligible compared to the kinetic term, and the theory is 
almost free: quantum fluctuations are "frozen". As / decreases, quantum fluctuations start to 
dress the system, which will eventually describes the proper quantum theory when / reaches 
its physical value. This argument will be used to define the boundary condition necessary to 
solve the self-consistent differential equation that we now derive. 
From now on, we represent a derivative with respect to / with a dot. We have then 




(12) 




(13) 



where the latter derivative can be obtained from the bar action ([21), such that 




(14) 
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In what follows, we omit the subscript c/ for the classical field. From eq.dH]), we obtain then 

r^//.*..'..<^.,.-.)A|.|,,,„_i|_}. (15) 

and, using eq.(l9]), we finally find that the exact evolution equation of F with / is 

t^fj,^. om^e + / / ,^.^y - .)^^ (^) " . (16) 

where the trace contains all the quantum corrections, and is regularized by the cut off A. We 
stress that this equation is exact: it is a self-consistent equation for F, which is independent of 
any perturbative expansion. 

Next step is to consider an ansatz for F, and we take here 



F = J d'^x (J^d^ed^e + u{e) 



(17) 



which allows an /-dependent IPI potential U{6). It is then sufficient to consider a constant 
axion configuration: the evolution of the effective potential is then given by F = VU{6), where 
V is the volume of space time. We have thus 



d d f 5^F \ ^ f d^p d'^q —p^q^ 



-(5(^)(P + 9)e"'^^"^'^^ (18) 



dx^" dy^XSejOy) J (27r)4 (277)4 /V + 
such that, taking into account V = 5^^\p = 0), 

<i4rf,i<.)(,^,)_|_^(_g_) =v/|^j5-j|l^. (19) 

where a prime denotes a derivative with respect to the axion. The integration over p leads then 
to the following non-perturbative and self-consistent equation for the IPI potential 



u{e) 



16772/ 



(20) 



where A is the UV cut off used to regularize the theory. There are several solutions to the non- 
linear partial differential equation ( 120|) . but, as expected from the spinodal instability discussed 
above, a flat potential U^ac is indeed a solution, and satisfies 

A4 

U^ac = (21) 

We consider some key points at this stage: 

• The flattening of the potential cannot be found with a perturbative expansion, but only 
with a non-perturbative approach: quantum fluctuations are so dominant in such a sys- 
tem, that they wipe off the classical non-convex features of the bare potential; 



5 



• The cosines in the bare potential do not appear anywhere in the evolution equation fl2U]) . 
because the latter is a self-consistent equation for the dressed potential. The cosines of 
the bare potential will play a role as a boundary condition used to solve this differential 
equation; 

• The axion, originally a Goldstone boson, and acquiring mass after integration of QCD 
degrees of freedom, eventually indeed becomes massless, after quantization of the model 

(ED; 

• The flat potential U^ac is not related to some vacuum expectation value of the axion field, 
and therefore does not need to satisfy any specific transformation under parity; 

• The equation fl?I]) can be written as a time-energy uncertainty relation, /r = 1, where 

r = 327t'^. (22) 

Hence, although we are dealing with a field theory at equilibrium, one might interpret r 
as the time that quantum fluctuations need to flatten the potential. This interpretation 
gives a physical meaning to the quantity Umc/-^^- 



Boundary condition. The integration of eq. (pT!) necessitates a boundary condition / = A ^ 
U = U\, such that 

t,„. = f7, + |lj (23) 

In order to determine Ua, we use the following argument. When / = A, the bare potential 
expressed in terms of the canonically normalized field (p = A6 is 



n=l 

where 



v{<p) = «n(i - (cos er) = ^^<P' + o (^1 j 



(24) 



4 

m 



^nan- (25) 

n=l 

Assuming m << A, the interactions are negligible compared to the kinetic and mass energies: 
the corresponding dressed flat potential Ua is then given by its one-loop expression, which is 
exact if the bare potential fl2^ is considered quadratic. The flattening of the potential is a tree 
level effect, and therefore of zeroth order in h: without further corrections, this would lead to a 
flat potential interpolating the minima of the classical potential, at zero energy. The one-loop 
correction, of first order in h, is then identified with Ua, 

Ua = ^TT{ln[{-d' + V")6^'\x-y)]} (26) 

\ , / A^ \ 1 



A^ 






[(- 


647r2 




A^ 


/A 


_ n 


167r^ 


V m 



AS ; ^ 2 ^ A4 
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such that the flat potential is finally 

327r2 V m 



Uvac^^\n{^). (27) 



To be consistent with the interpretation in terms of vacuum energy, we need Uyac > 0, such 
that the following constraint should be satisfied 

/A > m^ (28) 

which is a property common to every axion model. Note that the naive estimate of the vacuum 
energy density p^ac, based on the bare cosine potential, is 

1 f d^k /~ ^ ^ ^ 

^^^=27(2^^ + 7^ "16^- ^^'^ 

Although the ratio fK/m? is in principle large, the result fl27|) is actually of the same order of 
magnitude than the estimate fl2^ . because of the logarithm of this ratio. However, for heavy 
axions [13], the ratio fK/m? can actually be of order 1 given m = m{f), see for example [H] 
where m is quadratic in /. In such a situation, the vacuum energy ( 127|) can be much smaller 
than the naive estimate ( 129|) . such that the contribution of the axion vacuum energy is not 
necessarily orders of magnitude larger than observed and estimated values for Dark Energy. 
We note attempts to link the axion with Dark Energy [TJ] and that our result (I27p can be a 
model-independent starting point for such links. 

Interactions with other particles and discussion. One could argue that, because of the 
flatness of the potential, the initial motivation for the axion is lost, since no restoring force 
traps the axion in some vacuum expectation value, which compensates CP violation. In order 
to prevent this, we now consider a more realistic scenario where the axion is coupled to leptons 
with mass /i, via the term [16] 

gfoWi^ (30) 

where g is a dimensionless coupling, and we will discuss the conditions under which the spinodal 
instability can be avoided. Integrating out leptons, we obtain the effective potential (before 
quantizing of the axion): 

VefAe) = X: a.(l - (COS ^) ) - 2 y 1- [~~fT? ) ' ^''^ 

where the integral is regulated by the fermion cut off M. The spinodal instability, present with 
the cosine potential only, is avoided if the second derivative of the effective potential ( 13T|) is 
positive, such that the axion does not need to become massless. The original cosine potential is 
concave at = vr for example, where we evaluate the second derivative of the effective potential 
(ED): 

TT = —m — 



de^ ' ' 167r2 M2 + + ^2^27^2 

Q'P, 2 , o 2x2 2^, (. , M 



+ f^(/^^ + 3.VV)ln^l + -,^-^) (32) 
and we are left with the following situations (we neglect the mass /i): 
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a M » g f , the second derivative is 

Vi'n - - ^ < 0. (33) 



and the spinodal instability is not avoided. This is expected for a large fermion cut off, 
since fermions then have a negative contribution to the axion effective potential, and do 
not improve convexity; 



K'' ^ -m^ + cV/^ (34) 



if M ~ gfT^i then 

where c = [(3 In 2 — 2)/16]^/'^ ~ 1/4, and the spinodal instability can thus be avoided 
if gf is larger than the threshold ~ 4m, which is realistic, given the phenomenological 
constraints 



Note that, if the axion interacts with a scalar degree of freedom $, via the coupling //^6'^$^/4, 
where / is a dimensionless coupling, the correction f l5^ becomes 

and can be positive, if / is not too small, in which case no concave feature needs to be compen- 
sated by quantum fluctuations. This discussion can be extended to other interactions with the 
Standard Model, for example by coupling the axion to the photon via the interaction 

^OF'^'-FP'^e.^p,, (36) 

where h is a. dimensionless coupling. 

The essential results from this study can be summarised as follows. 

• When one takes into account quantum fluctuations of the axion, spinodal instabilities 
result in a flattening of the effective potential and a massless axion. We derive a model- 
independent form for the resulting vacuum energy which need not be orders of magnitude 
higher than estimates for dark energy. 

• In this quantized axion model, when axion couplings to other particles are considered, we 
find that spinodal instabilities can be avoided, allowing the axion to acquire mass as is 
the case with conventional axion models. 
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